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ABSTRACT 

Formulas are obtained for the mutual inductance of two helical wires whose 
axes are parallel. The formulas are given as series in oblate spheroidal harmonics. 
The principal part is the inductance of the two corresponding current sheets, the 
correction terms being small when the coils are closely wound, so that the pitch 
is small compared to the diameter of the cylindrical form in each case. 
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I. INTRODUCTION 

The principal part of the mutual inductance between two heKces 
whose axes are parallel is that between the two associated current 
sheets. This approximation may not be sufficient when the inductor 
is used for absolute electrical measurements of high precision. In 
this paper formulas are obtained which, in addition to the principal 
term, contain correction terms of the first and second order in the 
small quantities, such as the ratio of the pitch of the windings to their 
cylindrical radii. 

The electromagnetic force between currents in the helices may be 
obtained by differentiation with respect to their relative coordinates. 
Such helices are not used as often as coaxial ones, although they do 
occur, for example, in the Ayrton-Jones type of current balance as 
used by the English National Phj^sical Laboratory. 

II. DEFINITIONS AND INTEGRAL FORMULAS FOR THE 
MUTUAL INDUCTANCE OF TWO HELICES AND OF 
THEIR CURRENT SHEETS 

The two helices are shown in figure 1. The axis of the first, hi, is the 
a:-axis; that of the second, /12, is the line in the xy plane parallel to the 
X-axis at a distance r from it, where /•>0. If r=0, the helices are 
coaxial. 

619 
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The cylindrical radius of the first is ai ; its pitch, 27rpi ; and its axial 
length, li=2TrjpiNi, where A^i, the ''total number of turns," is not 
restricted to be an integer. The ''number of turns per unit axia, 
length" is 7ii = l/27r^i, not necessarily an integer. The initial plane! 







x=Xi^, contains its initial point, P^, and the end plane, x=Xe^ its end 
point, F^i, so that li=2TrpiNi=Xe^—Xi{^0, The current sheet, Si, 
"associated with the helix, /ii," is the cylindrical surface (generated by 
rotating the helix about its axis through an angle 27r), upon which 
there is a surface current of unit strength (in electromagnetic cgs 
units), uniformly distributed, the direction of the current at any point 
being that of the helical filament passing through that point. This is 
equivalent to defining the current sheet as that part, included between 
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the terminal planes x==^Xi and x=Xej of an endless cylinder which is 
completely covered by a helical tape winding of axial breadth equal 
to the pitch, so that no insulating space is left between successive 
windings. The angular current is Ui— 112^^1 per unit axial length. 
The axial current is l/27rai per unit length along a circular arc for 
which X is constant. 

If subscripts 2 refers to the second helix, h2, or its sheet S2, it follows 
from this definition that the mutual inductance. Mo, between the two 
current sheets is given by the double surface integral taken over the 
two sheets. 

where 

and R12 is the distance from a point Pi(xi,yij2i) on Si to P2fe,2/2>^2) 
on ^2, so that 

R2(x2—Xijdi,e2) = (x2—Xiy+al+al—2aia2 cos (do — di) 

+r^+2/'(a2 cos (?2— ^i cos ^i).] (3) 

The name ^^ mutual inductance^' is used here for the Neumann integral, 
so that when it refers to an unclosed circuit the implication is that the 
mutual inductance of the part (lead wires, etc.) which must be added 
to close the circuit must be evaluated by a Neumann's integral. 
With this understanding, it is evident that the term Pi^2Mi^2 in the 
numerator of eq 2 gives rise in eq 1 to the mutual inductance of the 
axial currents in the two sheets, whereas the term '^cos (62—61)^^ gives 
rise to that due to their angular currents only. 

In the integral (eq 1) the two independent variables, Xi, 61, are 
cylindrical coordinates referred to the axes of the first cylinder, and 
X2, 62 are cylindrical coordinates referred to the axes of the second 
cylinder, so that the angles ^1 and 62 each range over the interval from 
zero to 27r. If Xi, yi, Zi are the rectangular coordinates of any point 
Pi on Si, the position of Pi may be specified by assigning the pair of 
cylindrical coordinates Xi and ^1 (with cylindrical radius ai), or it may 
be specified by assigning the pair of semihelical coordinates Oi^, 61, the 
equations of transformation being 



\xi = Xi,+pi(di-di)] 0^9u<:2t 

\yi=ai cos di h where L ^. ^ 1_. a^o^m 

\zi=ai sin di J \d,^^di^de^=d,^+27rNi 



(4) 



Considering the pair of variables 61 and 6i^ for this range as two 
independent parameters which suffice to determine the position of 
Pi on Si, eq 4 are the parametric equations of the cylindrical surface 
Si, The locus of the equation ^i^— constant, is a helix, hi, of pitch 
2irpi lying on Si, wliile the locus of the equation ^1= constant is a 
straight line of length 2tPi (or less) lying on Si. 
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When the point Fi{xiyiZi) is also on the hehx hi, this may be desig- 
nated by Pni (Xh^j y^u 2hi)j the x^^, Vhv and Zn^ being given as in eq 4 
with the understanding that ^^ is a given constant and not a variable 
parameter. 

The equations 



X2=Xi^+'P2{B2—BiM 
y2=r-{~a2 cos 62 [, where- 
Z2= ^2 sin ^2 



(5) 

ei,^d2^de,^di, + 2TN2\ 



are the parametric equations of the sheet 5^2 when 62 and dt^ are both 
variables. The locus of the equation ^^g— ^^^^tant is a helix A2 of 

pitch 2x^2 lying on S2. 

The mutual inductance (Neumann's integral) of the helices, hi and 
^2, is 



• \ dsi \ d< 

Jo Jo 



Bv 



taken along both helices, where Si and ^2 denote total length of arc. 

From eq 4 it is found that the element of length dsi has the magni- 
tude 

ds,=^{al+fl)de, (6) 

and the direction cosines 

dxi_ pi dyi__ ai sin 61 dzi_ ai cos 61 ,^. 

The total length of the helix is Si, where 

5i=27riViVS!+F?=|;V^FFF?^ (8) 

where Zi is its "axial length. '' 

The same formulas, with change of subscript, refer to A2, so that 

dsijds2 cos (dsijds2)=aia2\ cos (^2—^1)+'^^ \dd1dd2' 

L ttia2j 

Hence the mutual inductance of the two helices is 

^aia^i'^deA ''de2F{xn^-XnJifi2), (9) 

where x^^ is a function of ^i only and Xn^ of ^2 only, given by eq 4 and 
5, the di^ and Bi^ being constants. 

In the general case of helices whose axes are merely parallel, not 
necessarily coaxial, and whose termina-l planes are arbitrary, eq 1 for 
the current sheets is sufl&ciently complicated, but eq 9 for the helices 
is even more troublesome. So far as we are aware, neither of these 
integrals may be evaluated in finite terms in general, although the 
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former may (in the case of coaxial current sheets) be so evaluated in 
terms of the three complete elHptic integrals. This is Jones' formula.^ 

The present paper is limited to what may, for brevity, be called the 
practical case, in which the two current sheets do not cut or touch each 
other, in fact they do not lie at any points close to each other. More 
precisely, considering the pitches 2'K'Px and 2Tr'p2 as infinitesimals of 
the first order compared to finite lengths like ai, ^2, then the distance, 
i?i2 of a point Pi on Si from a point P2 on S2, is, for all points, Pi and 
P2, finite (not a small quantity of the order of the pitch) . The various 
expansions in infinite series to be obtained below for the current sheet 
are not subject to this limitation but cover the ^^ impractical^ ' cases 
also. The restriction is necessary, however, in order to obtain a 
practicable quantitative expression for the relatively small difference 
between the mutual inductance of the two helices and that of their 
associated current sheets. 

lif(x,6) is a continuous and single-valued function of position on a 
cylindrical sheet S, it must be a periodic function of 6 with period 27r, 
and the line average of / along a helix on S is obviously very nearly 
the same as its surface average over S, provided that the pitch is 
small. 

In the following section the two are compared, and an expression for 
their difference is obtained to the second order of small quantities 
inclusive, and application is made to the problem in hand, because, 
with closely wound coils, the precision of electrical measurements at 
present attainable requires formulas valid to the second order. 

The general case of coaxial helices has been treated in an earlier 
paper,^ so the emphasis here is upon the noncoaxial case. 

III. DIFFERENCE BETWEEN THE FORMULA FOR THE 
MUTUAL INDUCTANCE OF TWO HELICES AND THAT 
OF THEIR CURRENT SHEETS 

Certain problems relating to the helix could be treated more ele- 
gantly by the use of semihelical coordinates or by orthogonal helical 
ones. A plane rectangle ruled with two mutually orthogonal families 
of straight lines could be wrapped around a cylindrical surface, and the 
parameters defining these families would be orthogonal helical coordi- 
nates. However, it has been found by trial that the comparison of a 
line integral along a helix with a surface integral over its sheet may be 
made most simply by the use of ordinary cylindrical coordinates 
(x, 6) because the x limits and 6 limits of the surface integral are inde- 
pendent of each other. 

As the prehminary formulation applies to either helix hi or ^2, the 
subscripts may be omitted. In figure 2, with x as abscissa and aS as 
ordinate, the rectangle represents the development of the cylindrical 
sheet S. The helix h is the system of parallel lines whose separation 
in the x-direction is Ax = 27r^— the pitch of the hehx. The slope of 
these straight lines in figure 2 is p/a. 

The variable 6 lies between zero and 27r, and the initial azimuth, 9i, 
also lies in that interval, but the end azimuth, de, as defined in eq 4 or 
5, does not, since de=di-\-2TN. We, therefore, (for the moment) dis- 

J J. V. Jones, Proc. Roy. See. (London) [A] 63, 198 (1898). 

3 Chester Snow, Mutual inductance and force between two coaxial helical wires, J. Research NBS 32, 26ft 
(1939) RP 1178. 
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tinguish between de and 6^,^ where ^^=9,— 27riV'=:^i+27r {N—N'), 
where N' is the largest positive integer which is less than or equal to 
^e/27r. When Be is thus chosen, it lies in the interval, 0^5e^27r, as 
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indicated in figure 2. The points Ai, A2, . 
Xn, where 

Xn=Xi—pdi+nAx for n=l,2j3j . 
Hence the equation of the helix for figure 2 is 
- =Xi+p(d—9i) when Xi^x<^Xi and ^,:<9<27r 

=Xn+pd when a:„<x<X;^+l and 0:<^<27r for 7^=1,2,3,. 

=Xe+p{d—ee) when XNr^x<^Xe and 0^(9<5c. 



An' have abscissa 



Xfi 



.N'-l 
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IjQii j{xfi) be a finite, continuous, and single-valued function of posi- 
tion on S. It is, therefore a periodic function of its second argument 
6, with period 27r. The same character is assumed for f\x,d) and 
j'\xfi), where primes denote partial derivatives with respect to the 
a:- argument. 

Hence, by eq 6, the line integral along a helix oi J{xfi){a^-\-p'^)~^ is 



V a2i-2?2jo je^ 

jdi Jo n=l 






+ jJiXe + 'p[d-de],d)dO^ 

The terms in the summation represent the jjrincipal part, the first 
and last integrals being relatively small (of the first order). Hence 
in them the integrand may be expanded by Taylor's theorem as a 
function of its a;-argument, about the terminal values of Xj so that 

-r^, rKx,e)ds=^a rdd'^'f(:xn+pd,d) 

■y/a^-i-'P^JO Jo n=l 

+a£dd[f{xe,e) +p(d-d,)f (xj)] 

+a\ dd[f{x„d)+2){d-e.^f'(:x,,0)]. (lO) 



To transform the summation, consider the following integral taken 
with respect to x, holding d constant, 

rxs' N'-l rxni-Ax 

fix,e)dx-=T, fix,e)dx, (11) 

J xi n^l JXn 

where the integral is taken along the horizontal dotted line passing 
through Ph and P in figure 2, the latter point P being the general 
point of integration with abscissa x=Xfy-\-{x—Xn), so that expanding 
f as a function of x (with 6 constant) about the point P/, on the helix 
in this n^ strip, we find 

J'*Xn-{-Ax f*x„+Ax 

j{xfi)dx= j{xn + [x—Xn]fi)dx 
Xn JXn 

= ^Jx'^\j{Xn,e) + {X-Xn)f ixnfi) +^^^/'(x„^)]. 

On performing the integration with respect to x (^-and x^ being 
constant), these become, after placing Xn=^Xn-\-'pd, 

rxn+Ax 
\f{x,d)dx= 

Jxn 
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so that eq 1 1 becomes 

rxN' 
f(x,d)dx= 

Jxx 

jv'-l 

[/(x„+2)e,e)+i?(^-e)/'fe+i7e,9)+|'[(,r-e)^+|^]/'(x„+pe,e)j- (12) 

Applying the same formula to i'ixfi) gives, after multiplying by 

^2Trp'^\pi7^-e)f^x„+pd,e)+p\'^-e)Y{x„+pd,e)], (12') 



which carries the approximation to the same order as in eq 12. Also, 
by eq 12, we may ^vrite to the same approximation 

p' (y'{x,e)dx=2wp''-z:pJixn+p6,e). (12") 

Jxi 71 = 1 

Subtracting eq 12' from eq 12, and making use of eq 12", gives 
On integrating the second and third terms, this becomes 

+£[(^-^^— 3-][/'fe',^)-/'fe,^)]. (13) 

In the terms of smallest order we may place f(xif'jd)=f(xe,d) and 
f(xi,d)=f(Xiyd). In the terms involving /(xj^/,^) and f(xi^d), which 
relative to the principal term, are of first order, we may place 

fiXN'yS) =f{Xe-pde>0) =f(Xe,e) —pdef {Xefi) 

f(x,,d) =f{Xi+p[27r^d,],d) =f{Xi,e) +p(27r-e,)f (x,fi). 
Also 

f(x,d)dx=^ I "^"^ \f(Xi.^) + ix-Xi)f{xM 

Xi Jxi I J 

= -(^^)[/(^*'«) +1 (27r-eO/'(x„fl)]- 
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j{xJ)\-{x--Xe)f{xJ)\ 
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Similarly, 

~— I J(x,6)dx=~ — I dx ' 

=-£[/fe,«)-%'fe,e)]- 

By use of these expansions, oq 13 becomes 

S/fc +pe,e) = ^ f%,B)dx 

The series on the right side of eq 10 is found by multiplying this by 
add and integrating from ^==0 to d=2T, Equation 10 then becomes 

To reduce eq 14 to a more convenient form, consider the function 
0{di—d), which is an odd periodic function of its argument with period 
27r, having a finite discontinuity when its argument is zero or 27rn. 
It is defined for all other values of its argument by the series 



0(6,-6)==^ 



sin n(di—d) 



71-1 n 



so that 



e(^,-^) = |-(\^), when 0<^,-^<27r 



(15a) 



(15b) 
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Applying the theorem (eq 16) to eq 17 gives 

-2a,^'jde2[v2e{e,,-e2)f{x,j2) +\vie'{o^-e,) -^]/(x,2,^2))- (i9) 

We may now introduce into the various integrals of eq 19 the ex- 
pressions iorf(x2,e2),f(Xe^,62)yJ(^i2^^2),f(Xe^A), and /' (3:^2,^2) obtained 
from eq 18. We retain only three orders of magnitude. 

Since/' occurs only in the second-order terms of eq 19, we need 
only the principal part of f ', which is the X2- derivative of the principal 
term in eq 18 so that we place 



'(X2A)=^\ dd,\ dxiF'(x2~XiAA) 

__^ f27r r 
'"2^Jo ^^ij^(^2-:^.p^i,^2)-i^X2^2-a:,i,6/i,0o) 



This is required only for the constant values X2=Xe2 and X2=^Xi^. On 
making these substitutions and taking account of eq 1, it is found that 



^=^o+T-ibrl"f^i2+^^2i+ 



47rViP2[ 

-V(j^2{Xe2iifiilfie^—Oi2{Xe2eifieifie^ 

+ Oi2{Xi^ev^eifii^~i^2{Xt^iJiV^i^ }, (20) 

where 

2a,a2^'dd,^dd2^2e{et-e^^ -p,Q {St, - ^i) J - '"'^^l'^ ^'^ ^ixt^tJiA) (21) 

and 

xtiti=xt,—xtu (22) 

where the subscript t suggests ^'terminal/' so that Xt2n is any one of 
the four possible distances between terminal planes, one of which is 
a terminal plane for the first sheet, the other for the second sheet. 

Each of the four x-differences represented by Xt^t^^ namely, 

Xe^i^^^^^Xc^ Xi^y Xe^e^^^^Xe^ Xg^j Xf^e^^^Xi^ Xg^, Xi^i^=^Xi^ Xi^ \^ZZ ) 

may be positive or negative. The notation dt^ indicates which ter- 
minal azimuth of No. 1, and Ot^ which terminal azimuth of No. 2 is 
associated with a specified x-difference. 
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The second-order term, 0^2, is an even function of its ar-argument, 
but the other terms are not so simple for 

J Jsi Jo 

[P2e{ee-e2)Fixe-xuei,e2) -P2O(et-e2)FiZi-xi,0iA)] (23a) 

and 

Wn^ ( fdS,2a, (ddi 

J Jsi Jo 

lvie{ee-e0Fix2-Xe^AM-pS{di-0i)F{x2-Xi^AA)l (236) 

If the integrals were left in this form, it would be necessary to dis- 
criminate between the various cases arising out of the possible com- 
binations of positive and negative signs among the four a:-values. 
The formulation of these cases would be tedious and awkward because 
of their lack of symmetry. 

This may be avoided without any loss in generality by transforming 
all the ^^-integrations of eq 23 into others in which the a:-argument of 
F has the same algebraic sign throughout the range of integration. 
To do this, we may write eq 23a in the form 

ddi I dd2\v20{de-e2) I dx,F(:x,—XifiM 

f" 1 

The first a^i integral is —•V2^{^i—^2) I dxiF{Xi—Xifiifi^\ 

I F{xe-xi)dx,= I F{x')dx'- I F{x')dx\ 

J^h Jo Jo 



The second is 



so that 

Wu=2aia2 



F(x,~xi)dx,^ I F{x')dx'- I F{x')d,x\ 

^h Jo Jo 

do, I cf6>2t2e(^e -^2) I F{x',dMdx'~ I Fix'fiMdx' 
Jo [ Lt/0 Jo J 

+:P2e(^i-^2) I F{x\diA)dx'- j Fix'^e.^dx' . 
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Similarly, 

2aia2 Cde, rS^2[2^ie(^c,-^i)[ f F (x', 6^,62) dx' - j^ F{x'AA)dx'j 

+Pie(ei-dJfjF(x\dMd^^ 

Adding these gives 

W,2 + W2l^ 
==(^l{Xe^i,,0i^,ee)-COi{Xe^e,A,A)+^l(Xi,^^^^^^ (24) 

where 

C^l(Xt,t^yOt^,et) = 

de, I de2[p2^{0t~e2)-vM^r~e,)]\ F(x',eMdx\ (25) 

This definition makes coi an odd function of its x argument. 

The mutual inductance of the current sheets ma^^ be put in a similar 
form by writing the a:ia>2 integrals of eq 1 

I dx, I dx2F{x2-x,)= I dx, I F{x')dx'-~ dxi I F{x/)d{x'). 

Xi, Jxij Jxi, Jo Jxi, Jo 

In the first of these double integrals let x''=Xe^—Xij and in the second 
let x^^=Xi^—xi. The first becomes 

Jf^Xei nXe^-Xl /•j.j-J.-j /*X' 

dx, I F(x/)dx'= I dx' I F(x/)dx' 
Xii J{) Jx'^'i Jo 

Jf^x,,ii nx" nx.,., r*x' 

dx" I Fixyix'- I dx' j F{x')dx\ 
Jo t/0 ,/o 

The second becomes 

J(*xe, nxi^-xi nxi^,, nx' r*xi^i, nx' 

dxi j F{x')dx'= I dx' I F{x')dx'- dx' F(:x')dx\ 
^•1 c/o i/o Jq Jo Jo 

Hence eq 1 gives as the mutual inductance, Mo, of the two current 
sheets 

Mo = ^oJ cOo(3:e2/l) -^0(3^02^1) +<^q(3^Vi) ~^Q(^<2^l)f ^^^^ 

coo(a;) =aia2 dd^ dB^ dx" dx'F{x'fiM, (26') 

Jo Jo Jo Jo 



where 



which is an even function of x vanishing with x. Its x derivative, 
Wo(aj), is an odd function of x vanishing with it. 
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An interpretation of the function, coo(x), may be obtained by con- 
sidering the special case where Xi^=Xi^ and Xe^=Xe2j which is 

where 

and ni = ll27rpi, 712=1/2x^2— number of turns per unit axial length. 
This shows that o)o(x) is one-half the mutual inductance of the two 
current sheets, when they have the same end planes, their common 
length being |.t|, and when each sheet has one turn per unit axial length 
(both axial and angular current in the sheets). 

Collecting these results, the mutual inductance, M, between the 
two helices (their Neumann's integral) is 

-o^ixi^.^^Oi^.e^^)}, (27) 

where 

co(x,2«i,<9«i,6',2)--coo(.T,2fi)+coi(a;,2^j^,^ii,(9,2) + co2(a;,2«j,^,^,(9,2), (28) 

where wo is given by eq 26\ coi by eq 25, and C02 by eq 21. The a:-differ- 
ences, x^g^i? between terminal planes of the current sheets and the 
terminal azimuths, 6t^ and 6t2, associated with each of these are given 
in eq 22 and 22^ It is now apparent that eq'28 is general, whatever 
combination of positive and negative signs the four x-differences may 
have. These signs are automatically taken care of in the definitions, 
eq 26\ 25, and 21, of cco, coi, and C02, respectively, which shows that 
and cuo and C02 are continuous even functions of x, wi being an odd 
function of x. 

It is evident also from these definitions, that coi and W2 (and hence w) 
are periodic functions of the terminal azimuths 6t^ and dt.^, since 6 is 
periodic with period 2tt. The period of oj will be seen to be 27r for each 
angle. Consequently, c^{xt^t2f ^t^j ^t^ is developable in a double 
Fourier's series of sines and cosines of multiples of these two angles, 
the development holding for all values of the angles and the coefficients 
being functions of Xt^t^ (and of r, ai, a2) 

Since wq is the finite term, both the terms (cos (^2-^1) and PiP2/aia2 
(constituting the numerator of the expression (eq 2) defining F) must 
be retained in evaluating coq. But since coi is in general an infinitesimal 
of first order and coo of second order, the term ^1^2/^10^2 is negligible in 
the definition of F, which is used for evaluating coi and 0^2 (since terms 
of third order have been neglected) . 

The system of subscripts in the designation of the three associated 
variables {xt^t^, 6t^, ^t) has now served its purpose of indicating which 
pair of terminal azimuths is associated with a given 3::-difference. 
From here on in the study of the co functions we adopt the simpler 
notation (x, 61, 62) for these three associated variables and place here 
for reference, see eq 22 and 22', 

x=xtt, di^dt d2=et^ ] 

[. (29) 

C^{xA,e2)=0^{Xt^t^,dt^,dt)\ 
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I This change necessitates a change in the notation of the variables of 

integration so that eq 26', 25, and 21 may now be written 

cooW =a.a. f'dx" n.' rde.' fW^^^^f^-f^ (30) 

Jo Jo Jo Jo -ni2ix ,C/i ,t/2 ) 

2aia2 Cdx^ rde,' rddAv2Q{d,-e,')-vS{e,-e,')]^^^^^^ (31) 

Jo Jo Jo Ui2\X ,Ui ,d2 ) 



+f\^^\Q2-e',)~-^--2^,V20{d,-d[)Qie-d'^ 



icos ((9;—^;) 






Equation 28 becomes 

In the relatively small terms wi and C02, a sufficient approximation to 
the exact expressions 

m-^'"^ and em-ll=±'^P (34) 

is 

G(^)=^(^sin ^+^ sin 2^) and e\e)~^^=-'^i^^os ^+^^) (35) 

Using this approximation, eq 31 becomes 

coi(:r,(9i,^2)=^G^i«'2 I dx' \ ""dd'A rfd^^i sin (^i— ^O+^sin 2(^; — ^1) 

-;,{sin (.;-..)+^sin 2(^;-^.)]f£|?£f . (36) 
and eq 32 becomes 

co2(a;,e,,92)=^^^ ft p^pjfcos (&;-0i)+|cos 2(0;-9x)l 

+l)f cos (0^-62) +1 cos 2(9^-52)1 

_2^=[sin(e;-e,)+|sm2(e;-^,)] 

[sin(.-.2)+isin2(.-.2]}-^|yf (37) 
For reference is placed here an alternative approximation 
Q\e) = '^(\-sm |\ when 0<9<27r 

= ^n+sin^Y when — 2;r<(9<0 



(38) 



269047—41 2 
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which is seen to be a good approximation, since 

32/7r3- 1.032 



and 



.2 32 « sin (271-1)2 



Q'W=y ^S^.—7o^:z^^\^^ 1^ when 0<KW2. 

For the type of expansions required here the three-term approxima- 
tion for 6^ given in eq 35 will be used. It has about the same precision 
as the two-term approximation, eq 38. The latter is equivalent to 
one used ^ in considering coaxial helices, as it led to a simple expression 
for the term corresponding in that case to C02, which was in finite terms. 

A comparison of the approximate expression adopted for 9 id) and 
for &^(B) —IP/12 in (eq 35) with the exact series (eq 34), for these func- 
tions shows that we have retained only the first two terms of the series, 
and have arbitrarily inserted the factor 7rV9, which is slightly greater 
than 1, and in fact it is sufficient to replace this by 1 in the second- 
order terms. 

In the following section the w-functions are expanded in a double 
Fourier^s series in the terminal azimuth angles Oi and ^2. 

IV. EXPANSION OF CO IN A DOUBLE FOURIER^S SERIES OF 
THE AZIMUTHAL ANGLES 

In figure 3 is shown the projection of the two parallel cylinders 
upon the yz plane, and ri2 is the projection of the distance Ru from Pi 
to P2, so that 

R?2=^^'\-f'i2j where a: =0:2—2:1. (39) 

The X axis being taken as the axis of cyhndrical coordinates, ai and 
61 are cylindrical coordinates of Pi, but the cylindrical coordinates of 
P2 may be taken as rj and ^2, as shown in figure 3. Since ri2 is one 
side of the triangle P1O1P2, the other two sides, ai, and r2, including 
the angle ^1—^2, 

r?2^a\~2air2 cos ((9i-^;)+r^. (40) 

Since l/i?i2 is an even periodic function of the angle (^1—^2), its Fourier 
series may be written 

1 2 " 
w-=-j=T S €n«n(a:,ai,r2) cos n(i9i~(?;), (41) 

where €0=0 ^^^ «n=l, if ^t^O), and where 



and x, tti, r2, $2 are held constant in this integration. 



cc 



<l>n(x,aur2) _ 1 I cosn(di — e2)ddi 



(42a) 



8 Reference 2, p. 259-261. 
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Another integral representation of this function, which will be 
derived later and which involves BessePs function, Jn, is 



<t>n(x,aur2) _ 



V' 



aiTi 



= (p''^'Jn(ait)Mr2t)dt. (42b) 




FiGUBE 3. — Projection of the two helices upon a plane perpendicular to their axes. 
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This integral may be evaluated, the result being a function of the 
single dimensionless variable k, so that 

<i>n(x,a„T2) F"+i ^V+2) / 1 1 \ 

Vairo 27rV^i^2 r (271+1) \ 2 '2 ' > y v ^ 

where 

and where i^ is the hypergeometric function. It is evident that the 
function 4>n,(^—Xijrjri)/-y/rri is a solution of the partial differential 
equation 

(Dt+Dt+lD.-^^ = (45a) 

(dI^+D^,+Id.__-^):^=0. (45b) 

It is further evident that 0^ is a symmetrical function of the two sets, 
x,r, and Xi,ri, becoming infinite logarithmically when they approach 
coincidence. It is the Green's function and the nucleus of many 
integral equations which arise when a solution of eq 45 is required 
and which has assigned values (or assigned values for its ^'surface 
density of charge'^ on surfaces of revolution about the x-axis. Thus 
the function 

Vo{x,r) = (j>o(x—Xi,r,ri)/'ylrri 

is a potential at a point P whose cylindrical coordinates are x and r. 
It is symmetrical about the x-axis and is due to a unit circular line 
charge distributed with uniform linear density over the circular arc 
of radius ri in the plane Xi, which is coaxial with the a^-axis. 
More generally, if n is any integer, the function 

cl)n{x—Xi,r,ri) 



and also of 



V: 



is the solution of eq 45a at a point x,r, due to a unit circular line 
'^source" at Xijri, in the sense that the solution of eq 45a, whose 
surface density of source a is given on any surface of revolution 
(a:i=function of rO, is found by integrating 

r .(rO<^.(^-x.,r,rO ^^^,^^^^ 

where dsi is an element of the generating arc of the given surface of 
revolution, lying in a plane through the x-axis, and the integral is 
taken between the appropriate limits. If a is known, this gives the 
solution, but if the function is assigned, this is an integral equation 
to determine cr, after which this integral gives the solution. The most 
appropriate expansions of this nucleus of the integral equations in a 
given coordinate system are those which are the canonical bilinear 
expansion of the nucleus in terms of its normal functions for that 
system. 
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The role of this function </>„ for solutions of eq 45a, having assigned 
values on surfaces of revolution, is as fundamental as that of l/Rp^p 
in the general theory of electrostatics, and in fact its expansion in 
normal functions follows from that of 1/Rp^p by eq 42a. 

Some coordinates pertinent to the problem of this paper may here 
be mentioned. 

The first is the spherical coordinate system. If the values of a 
solution U of eq 45a are assigned on the sphere made by revolving the 
circle xl-\-rl=Rl about the a:-axis, and R, f, 6 are the spherical 
coordinates of P, while (i?i, fi, ^i) are those of Pi, where R=x^-\-r^j 
^=x/R and tan d=y/z (where 6 is the same as in cylindrical coordi- 
nates), then it is known that 

1 ta ° V(m TJ-J-I") 

= 2 Z:«. cos n{d-ed S w!^,:li P:,{OR-'"-'Pl(^r)Br (46) 



Rp^p ^ ' ^nT{m+n+l) 

when i?>i?i. (with interchange of R and Ri for the case i?<i?i). 
The expansion corresponding to eq 41 is 



so that 



RTI Jz=tl2^n(l>n{x-Xi,r,ri) cos n{d-e,), (47) 



(j)n{x—Xi,r, ri)__^ V{m—n hi) p^ />.Np-m-i pn /> \-pm /^o\ 



when jB>i?i (with interchange of R and Ri when i?<i?i), where 
-Pm($) ^ is the associated Legendre function (Ferrer's). This is the 
canonical bilinear expansion of the nucleus in terms of normal func- 
tions for the spherical surface^ i?i by means of which solutions of 
eq 45a may be found when their values or the values of their source 
density on the sphere are assigned. 

When the solution of eq 45a has assigned values on an oblate 
spheroid, oblate spheroidal coordinates are used. The bilinear expan- 
tion for this nucleus in such coordinates is used in a later section and 
applied to the limiting case where the spheroid becomes a circular disk. 

The expansion in toroidal coordinates is also intimately related to 
the present problem, the function (l)n(k) being a toroidal function 

(j)^(^Jc)=-Q^_J——l \ where Qw is Legendre's function of the second 

kind. 

The derivation of eq 42b may be made by use of Heine's integral 
formula, 

Using Neumann's addition formula, 

00 

=2J2^nJn{ait)Jnir2t) COS niOi-d'^) (50) 
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gives 

^=2X16. cos n{e,-e',) {h'^KJn{ait)Jn{r2f)dt, (51) 

xll2 n=0 JO 

which, by comparison with eq 41, gives eq 42b. In the expansion 
(eq 41) or its equivalent (eq 51), both r2 and $2 are functions of $2, 
Since r2 is one side of the triangle, P2O1O2 (fig. 3), the other sides, r 
and (7.2, including the angle {-^—62), 

rl=r^—2ra2 cos (tt— ^2)+g^2- (52) 

To obtain the Fourier's cosine series in ^1 and 62 for I/JS12, the term 
Jo {ri2t) in the integrand of eq 49 must be expressed explicitly in 
terms of these angles. Hence in the expansion (eq 50), the term 

Jni(^20 e'^^^^^ must be expanded as a Fourier series in 62, and for this a 
more general formula than Neumann's addition formula (eq 50) is 
required. 

Neumann's addition formula is a special case of a more general 
addition formula derivable from that of Sonin,* which, applied to the 
triangle P2O1O2, may be put in the form 

e/„,(r20e-*'*^^i- S Jn,-n,{rt)Jn,{a2t)e-'^^\ (53) 

Using this in eq 50 gives 

e/o(n20-= i: Z; Jn,-n,{rt)JnM^t)Jn,{a2t)x''-^\--^'^\ (54) 

Til = — 05 7?2= — CO 

By use of this expansion in eq 49, the double Fourier's series expansion 
of the function IjRn defined in eq 3 is found to be 

^ ^ S S <i>n,,nXx,r,aua2)e'^-^'^'--^'^' 



Bu(,X,Sif62) -yai^2 ni»-» nj=-oo ^' ^ 
1 



-7= S S (l>n.,no(Xyr,aua2) cos (:nA— 71262)9 (55) 



/(ll(Z2 711= -00 nx 

which may be written 



+S IZ 'i>n^n2,ni<iOB[(n-{-n2)ei-nM\, (56) 



or 
1 



B\2{xfi\fi2) 



4-) CD CP f 

~r=S Zle«i€J(<^non2 + <;^«it-n2) COS Tli^j COS 7l2l92 
VCtia2ni=0 nj=0 I 



+ (««i,n3-0ni,-«2) siu n^ siu Tljl^s • (56') 



« N. Nielsen, Haodibucli der Tbeorie der Cylinder Fun^tionen, p. 287, Teubner (1904), 
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The coefficients <^r?i, ^2 are even functions of x, which are represented 

by 

or by Fourier's rule, 

/ N Jchao HJa Pin COS (niOi—nodo) ,.-, v 

which shows that 0„j,_n2=</>-7ii.n2- Oil inserting the expansion (eq 56) 
m the three equations, 30, 31, 32, noting that 

Oiie—e'i) cos ne[dd[=- sin nOu if n^^O, =0 if /^--O 
9(^1— ^'i) sin ne[de[= cos ndi if tIt^^O, =-0 if n=0 

O^^i — ^'i) cosn^id^i=--2 cos 71^1 if n?^0, =^ ; if n=0, 

the Fourier's series for the function o) is obtained. The term coq 
corresponding to current sheets may be written 

coo=ajii+-^-cooo, (58) 

(iia2 

where the principal part of co is wn due to the angular component of 
current in the sheets. The term 

JhP2 

a.a,''^ 

is due to the axial components of current in the sheets and is an 
infinitesimal of second order. 

The functions o^n and coqo are special cases of the function co„j„2 
defined by 

^ni,n2(^/>^l>^2) =47rVfl^lC^2 I dx'' j dx' (i>n^,n2{^' ,r ^auO^ , (59) 

'ni,n2(a^/,c^i,tt2) =47ry^ j da:'<A«i,n2(aj',^,cti,^2) • (60) 

Jo 

The first-order infinitesimal coi given by eq 31 is found to be 

coi(x,^i,^2) == —Vi S — r-^ Wn+iA^yr,aiya2) +a)'n_i,-_i(a;,r,ai,a2)] 
71=1 n 

^Sin 7l(7r — ^2)r / / \ » / / M //>-i\ 

—V2 2Li [co'n+i,i(a:,r,a2,ai)+a)n-i,-i(a^,r,a2,ai)], (61) 



so that 

CO 



— 1 n 
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and the second-order term 0^2, derived from eq 32 becomes 

^02(^,^1,^2) 

==4.T^^|ala2\pl^^^^-^ 
[ n=l n2 

71=1 n 
00 00 "j 

+ PlP2^ S — --[(0.+l.-n,^-l + </);^^~l,-7^,-l) COS (7li0i + 7l2^2) 



— (*n,+l.n3+l + </>;^^-l,n,-l) COS (tIi^i — 713^2) ]j^ (62) 

where the arguments of all the </)-functions in the double series are 
(x, r, di, a2), respectively. In deriving these expressions from the 
expansion (eq 55), we have made use of the fact (shown by eq 57b) 
that (t)-n,-n —<t>n,n ^nd also of the relation (shown by eq 57a that 

<t>n^,n^(x,r,ai,a2)=(l)n^,nX,r,a2,ai) , (63) 

which obviously hold for the functions a;„^„^ and co\^ny 

The approximate expressions for 9 and 0^ given in eq 35 merely 
indicate that wi is found from eq 61 by retaining only the first two 
terms of each series and inserting the factor 7r79, while the second-order 
term, C02, is sufficiently represented by merely extending the ni, and 
712 summations each from 1 to 2 in the series given in eq 62. 

In the case of coaxial helices (r=0), eq 57a, shows that 

<l>n^n^{XyOjaiya2)='<l)n^ix,aiya2)dn^n^=(j)n^{k)8n^n^f 

where 

^^+(^1+0^2)^ 

and dn^,n='^y if 7ii=n2, = zero otherwise. Equation 61 shows that 
wi = 0, and eq 62 shows that 

0)2-^-4.'l^^PlP2^JalCi2^[c|>n-l{k)+c|>,,^l{k)] ^Q^ ^H 1^- 2), 

71 = 1 lO 

where 61 — 62 is the azimuthal difference of the two terminals. 

Hence the surviving terms are equivalent to those previously 
obtained for the coaxial case in eq 49 of the paper quoted (footnote 2). 
The terms 7ii7i2coii(x), nin2o^0Q(x), 711^2^2(^,^1 — ^2) of this paper correspond 
to o:o{x), o)x{x), Siiid o:a{x,6x){6x='-di — 02) of that paper, where 7ii7i2 

For the coaxial helices in a current balance, the term C02 compared 
with Wo was found to be 10 parts in a million. 

To derive a formula by which the mutual inductance of the two 
helices may be computed in the general case, where their axes are 
merely parallel, it is necessary to obtain expressions by which the 
functions (t)n n and their first and second partial integrals with respect 
to X may be computed. None of the integrals, eq 57a, 57b, 59, and 
60, may be integrated in finite terms, so we are driven to expansions 
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in series. Such expansions in oblate spheroidal harmonics are ob- 
tained in the following sections. 

By reference to the relation (eq 63) it is evident that ^^jj^^, where 
Til and % are any integers (positive, negative, or zero), is obtainable 
(by interchange of the arguments, ai and a2 when necessary) from one 
of the following three irreducible cases: 

Case Jl; 7^2=?^l>:0 
Case B; 0:<7i2<7ii 
Case C; 7i2<0 and 7ii>0 

However, this requires that we obtain expansions for 0„„ in these 
three cases, without committing ourselves to a restriction that a2<C^ij 
that is, one that is valid when ai and ^2 are interchanged. 

V. PARTIAL DIFFERENTIAL EQUATIONS AND BOUNDARY 
CONDITIONS SATISFIED BY THE FUNCTION (^„^,„^ 

From the differential equation satisfied by BesseFs function, com- 
bined with the integral representation of (t>nj,n^ in (eq 57a), it is found 
that 0„i,7?2 satisfies the two partial differential equations 



and 



\^Dl+Dl+lDr-^Jl^^=^0, where n==n^~-n,^i), (64) 

\Dl+Dl^+^Da~fi^'==0, (65) 

and also a similar equation, where ai and ni are replaced by a2 and %. 
As to boundary conditions, it is evident that .-^ vanishes when 

-yjx^+r^ -^oo , or when ^Jx^-\-al ->oo . The function is also regular in the 
neighborhood of r=0 or (ai = 0) and is finite, continuous, and single- 
valued in general and an even function of x. 

The partial derivative with respect to x, (t>n^,n^ is an odd function of x. 
Defining a by 

_ 4>l^,n,{+^,r,ai,a2) (/>„'^,„^(--0,r,ai,ao) 



*^7^'y/aia2 27r^aia2 

1 C"^ 

=2^ J^^^«,- n,(rO J«^ {a, {)Jn,{o.2t)dty (66) 

it will be seen that cr must be zero when r< | ai — % | and when r>ai + (i2, 
for the integral representation (eq 57b) exhibits ^n^jn^ as a double line 
integral around two parallel circles of radii di and ^2, x being the 
distance between their planes and r the distance between thei r axes. 
In the case ni = n2^l, the function </)ii multiplied by 47r^V<^i^2 is the 
mutual inductance of the two parallel circles. 



642 Journal of Research of the National Bureau of Standards [voi. t5 
Equation 57b may be written 



4:Tr^-\/aia2<t>n^,n^{x,r)= I aa cos 77.2^^2 I 
Jo Jo 



(Xi C OS nididdi 

~R^2 



J2ir /»27r 

a, sin uAdd, j ^^-^^cZ^i- (67) 

Each of these integrals is a special case of double-line integrals taken 
around the two circles of type 

'i>{x)^ ^Ue2)ds,^-^-^^§^, (68) 

where 

dsi=aiddi and ds2=a2dd2 

are elements of arc of the circles, R12 being the distance between the 
two points of integration, Pi in the first circle and P2 on the second. 
The first of radius ai may be considered in the plane x—0 coaxial 
with the .7:-axis, the second parallel to it has its center at a point whose 
cylindrical coordinates are (x;r,d). Since l/Ru is so fundamental in 
the theory of electrostatics, and because electrostatic problems are 
familiar to most of us, some guidance may be obtained by availing 
ourselves occasionally of an electrostatic interpretation of integrals 
involving the factor l/i?i2 in the integrand. Thus the first line integral 
in eq 68, which is 

J R\2 

is the electrostatic potential at any point P2 on the second circle when 
the first circle is a line charge with linear density fi(di). Hence, if 
the second circle is a line charge of linear density /g (^2), then 

-^'(a;)-.-~/?^0- \j2(02)E,ds2, where £^^--7>;F 

the a:-component of the electric field of the first circle at any point P2 
of the second. Consequently, —(l)^=Fx=th.e cc-component of the 
force which the first circle exerts upon the second. Clearly this 
vanishes w^hen x=0 if the circles do not touch or intersect in this 
coplanar position, that is, 

<t>'(0) = if 0<r<|ai~a2| (69) 

^0 if ai+a2<Cr<^co, 

In the remaining case, \ai—a2\<Cr<C(ii+(i2, the circles intersect if they 
are brought into the same plane while keeping r constant. This is 
the ca^e for which figure 3 is drawn, assuming a2<Cai (but the argu- 
ment may be made by reference to that figure, keeping in mind the 
other possibility, ai<a2). In this case the values of di and 62 corre- 
sponding to one point of intersection Pq are called 61 and 62, respec- 



f 
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Equation 14 may therefore be written 

-2aJJS^[^e(^,-^)/fe,^) + [i>202(^,-^) -^']/(^/,^) j} (16) 

where "Qe tas been replaced by Be, the distinction between the two be- 
ing now no longer material, since de appears only in {Se—d)^ which 
is a periodic function of its argument with period 27r. 

The parenthesis on the right of eq 16, which has the factor l/27r2?, 
contains three orders of magnitudes. The first, or finite term, is the 
surface integral of/ taken over the sheet S. 

If eq 16 be divided by 2TrNa, the first member becomes l/sj^'fds^ 
the line average of / taken along the helix. The principal term on 
the right becomes l/SJ^J^JdS, which is the surface average of / taken 
over S. Equation 16 is therefore the required relation giving the 
difference of these two averages as a relatively small quantity in which 
two orders of magnitude have been retained. 

By two successive applications of the transformation represented in 
eq 16 we may convert the two successive line integrals along the two 
helices (which represent their mutual inductance, M) into two suc- 
cessive surface integrals over their equivalent sheets plus small cor- 
rection terms. To begin this, eq 9 may be written 



M= 



^2 



-4? 
ijo 



dS2f(X2jd2) 



the integration being along the second helix, 



(17) 



f(X2A)=' 






F{x2~XiAM, 



the latter integration being taken along the first helix. 
By use of eq 16 

/fe,e2) =2^( JX dSiFix2-XiAA) 

+ 2ai[deJ^pM0e-Oi)F(x2-x,^,di,e2) 

-[pfe^(e.-eo-^']F'fe-x,,AA) 

-2ai(dei>^pM0i~9i)F(x2-Xi^,di,d,) 

-[_Pie'ie^-e,)~^']F'{x,-x,^,eM]). (is) 
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tively. The angle 2</' between the intersecting arcs may be taken 
as the angle O1P0O2, so that 

2yp^dl—el, where 0<^;<7r and 0^^^<7r, 1 

so that 0:<2i/'<7r ) ^ ^ 

The coordinates, 61=27:— &i and d'l=2'K—el refer to the other point 
of intersection Pq. 

When the two circles approach the same plane (x-^+0), the 
infinitesimal arcs in tlie neighborhood of Pq (or P'q) contribute a 
finite amount to their mutual repulsion F^ in the a:-direction, because 
R12 vanishes at Pq (and Pq). Arcs at a finite distance from these 
points contribute a vanishing amount to Fx. The curvatui-e of the 
intersecting arcs becomes irrelevant, but of course the angfe 2\l/ at 
which they intersect is important. The value of the integl'al may be 
computed as in the case of two endless straight-line charges approach- 
ing intersection, their linear charge densities being variable, but only 
the values of these densities at Pq (or P') being effective. If we 
allow for discontinuities in the functions /i (^i) and J2 (^2) (as for 
example in eq 31 and eq 32, in the integrals defining coi and C02), we 
find, when [ai — a2Kr<ai + a2 and (j>, defined by eq 68, that 

,// , nw -^ f/ /i(^M-Q)+/i(^;-0) V /2(^g+Q)+/2(^2"-0) \ 
|/ /i(^^?+Q)+/i(^^2"--Q) Y /2(^^g+0)+/2(^^^-0) \ 
|/ 7r-4^ ^r/ /,(^;+o)-/,(^?-o) y /2(^^fQ)--/2(^^~^ 

, / /i(^^;+Q)-/i(^^?-0) Y /2(^1+Q)-/2(^1-Q) \l ^^^) 
For the case where /i and J2 are continuous functions, this becomes 

-<^'i+^)-^2^U^^^^U^)-^:f^^^'^)Um. (71') 

Using this in eq 67 leads to the following expression for the o- defined 
in eq 66: 

(T^-O ifO<r<^2 

—.Tl^. '"^ y-^-^°) if A<r<M, (72) 



where 



n^{A\-Af) sin 2^p 

=-0 if^i<r<co 

L so that0<^2<^i. (73) 

■As^iai— aslJ 

The absolute value |ai— a2| is used for ^2, so that we need not commit 
ourselves on the question as to which of the two circular radii, ai or 
a2, is the larger. 
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The second member of eq 72 is a function of r, since 61 and 0% and 
2^=dl-~di are functions of r. However, it is more convenient to 
reverse this and to express a and r as functions of xf/, which may be 
done by expressing ^J as functions of \^. 



To do this we define |o by 



t — ^/i t?_2V5^__ /1 M2 



ai+a2 



0<&<1 



Considering the triangle O1P0O2 of figure 3, it is evident that 
r={al~-2aia2COs2\l/+al)^=Ai(l — ^l cos^ xp)^ 

=»<'-S«")'('-5;«"*)' 



sm 



^^^ZF::aiV(^?~r^)(^^-^i) 



cos C/i = 



%— a2 cos 2i/'] 



sm 6>J 



,0 ^2 sin 2i/'| 



and < 



cos (92 



_ai cos 2\l/~a2] 



sin eo^^!LJHi2^ 



From tliese relations it follows that 



e»?=- 



1 ^0„_2i^ 






.i^;_ «ie''''^-g2 . 



:e«(«?+2«. 



(74) 



(75) 
(75') 
(76) 

(77a) 
(77b) 



When considering di and 61 as functions of r, ai, and ^2, this may be 
indicated by writing ^?=^?(r,ai,a2) and 62=^62(r,aija2). It is evident 
from these definitions that when we interchange the arguments ai 
and 0^2, 61 goes into tt— ^2 and ^2 goes into tt—^?, that is, 



and 



^?(r,a2,ai) = 7r-6'^(r,ai,a2) 
(9^(r,a2,ai) =:7r— ^f(r,ai,a2). 



(78a) 
(78b) 



On the other hand, 2\l/{=6l—6i) is unaltered by interchange of argu- 
ments ai and ag. 
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By eq 77 we may write (if ^ ( ) denotes the real part) 



A 2 / ni-m 



cos (■/ii^;-n2^) = 6i^( 



/l+I^Y'-"' 



„,_-^r(n,-n,+i). 



cos 2(^+nj)^ 






Til — 7l2 

(l-^?cos^ -A)-^ 



(-»t:i)' 



COS 2A:i/' 



Hence, if we define the integer n by 

n=ni — n2j so that n^O in all cases, ^, B, and C, (79) 

while 7^2 nia}^ be positive or negative, then 



(7=0 when 0<r<^A2 

2 /i+ioV r(7i+i) 

|l7r2A2sin2;/'V 2 / ,, .^ 2 , « 

(1— JoCOS^\^)2 

cos 2 (A- +^2) 'A 



2^<-'><T+i)"i 



(80) 



Tik+l)r(l+n-k) 
when ^2<C^<C^i 
— Q when Ai<C^^ 00 

In case n=^0, this series reduces to the simple term cos 2n2\p, so tha 

_ 2 cos 2m\^ 
''"'^~x2(^f-.4i) sin2iA' 

Quite a different electrostatic interpretation of </>„ ^ will now be 
useful. If X, r, and ^ are considered as the cylindrical coordinates of 
a point P in space, the function V(x,r,d) = (aia2)~Hni,n2{^,'^,cii,oi'2) 
cos {ni—n2}Q is the Newtonian potential at any point P due to a simple 
distribution of charge upon the annular surface {A2<jr<^Ai,x=Qi), whose 
surface density is 0- cos {ni — n2)d, so that the a defined by eq 66 or 80 
plays the same part for solutions of eq 64 that surface density on the 
annulus plays for solutions of Laplace's equation. 

In the following section an expansion for <i>n^^n2 is obtained in terms 
of the normal functions of oblate spheroidal coordinates, these being 
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used because in terms of them 4>n^,n2 admits of a single expansion 
valid for all real values of x and for all positive or zero values of r, ai, 
and ^2, whereas ordinary spherical coordinates require three different 
expansions according as 0:<r:<x42, ^2^^=^^i, or ^i:<r:< 00, and the 
expansion in tiie second of these cases is somewhat complicated. In 
this respect, toroidal coordinates are superior to spherical but inferior 
to spheroidal, as the expansions require only two cases. 

VI. EXPANSION OF t^n^n^ IN OBLATE SPHEROIDAL 
COORDINATES 

Instead of the cylindrical coordinates {x, r, 6) the (oblate) spheroidal 
coordinates, ^, ??, 6 (the same d), me,y be used to designate a pomt in 
space, where 



x^Ai^rf 
Solving for J ad 77 gives 



, where 




^'^2Zfn"^'-^'+V[^^+(r+^40^][x^+(r-A)^] 



1 



'' 2Al{ 
Also from eq 81, we find 



r'-Al+^[x'+(r+A^yW+{r-A,y]\ 



(^xHrfr^==.4?(r+.'0[i^2+if^2]- 



(81) 



(82) 



(83) 



The locus of the equation, r; == constant, is the entire surface of the 
oblate spheroid, whose equation is 

x'/AW+r'/Al{l+r]'^) = l, 

whose foci are at (x=0, r= ±Ai). In particular, the locus of r;^0 is 
both sides of the circular disk (x=0, 0:^r, :^Ai), at the edge of which 
^=0, and ^ increases from zero to 1 on the side x—-\-0 as the point 
passes from the edge to the origin, or a:-axis, while ^ decreases from 
to —1 on the side x=—0 as the point moves down to the axis. 
The locus ^=^ + 1 is the positive x-axis, and ^= — 1 is the negative 
a;-axis. 

The hyperboloid which is confocal with every such spheroid whose 
equation is 

-xyAie+r'/AKl-e) = l 

must be considered as the two loci, that half for which x>0 being 
the locus ^=|i>0, and that for which x<CO being the locus ^=— fi<CO. 
Large values of rj correspond to points far away from the origin, since 
l-{-r}^~{x^-\-r^)/Al^-^'\ The differential equation (64) becomes 



D,[(i-e)D,]- 



n^ 



nj n2 



1 — $1-V^1<^2 



^^D,[(l + yj')D,]- 






(84) 
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which has solutions of the form 

^^PiORiv), (85) 

where P(£) is a sohition of 

A[(1-J2)P']+[m(m+1)-Y^2]p=0, where P^^D.P (86) 
and R{rj) is a solution of 

D,[(l + ,^)i?']-[M(M+l)-i^]«=0 (87) 

n=?ii— 7i2. (88) 

As shown in the preceding section, it is sufficient to consider ti— ni— 712, 
where n is a positive integer, or zero, while 712 is either positive, 
negative, or zero, and ni'^0. Equation 87 goes into the same type 
as eq 86 in the independent variable 2, where r}=iz. 
One solution of eq 86 for the case where — 1<^<1 may be taken as 
P=Pl(^), which is (Ferrer's) associated Legendre function. This 
reduces to the accepted value when j^ = an integer, if we define it, for 
general values of ijl,v when|l — J|<2 by 

iM(0"2T(M-^+i)r(H"iy V~''''+^+^''+^'"T"/ ^^^) 

We make use of the two recurrence relations 

(2M+i)eP:(r)=(M-^+i)n+i(f)+(M+^)p.-i(f) (90) 

(2;.+ l)(l-f^)P;'(?) = -(^-.+ l)MP,;i(e + (M + ^)(M+l)PM-i(i^-(91) 
Also when n is a positive integer 

p:a)=ii-e)^mp.{o. 02) 

We require also the following expressions when /*=v+s, where s is a 
positive integer or zero. 

r-.i^Ki) 2" r(s+l)r(v+l) 

F(^-s,s+2v+l,v+l-^y if |1-||<2 (93) 

P:..(^)- ^-^^'y-^^^^^^-+:+^> F i-s,s+.+i,m (94a) 
Vtt r(s+l) 
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Also when |1 — ^^Kl, there are the two relations: 

p. ... 2-(l-a^^ V{s + v+\)V{s + v+^) 

^^+2au; r(^_j_i) r(s+i)r(s+i) 

i^(-s,s+.+i,.+ l;l-a (95a) 
p. ... 2^{l-e)iv{s + v+l)T{s+v+i) 

^.+26-+iu;- Y{v+i) r(s+i) r(s+|) 

^i^(~v+^-fV+i;i-a. (95b) 

The hypergeometric functions in the last five expressions are polyno- 
mials, if s is a positive integer or zero, but eq 93, 95a, and 95b are 
valid if V and s are unrestricted. 

The importance of the case ix = v-\-s lies in the fact that if v is any 
given constant greater than —1, the infinite set of functions PIj^^^^ 
constitute a complete set of orthogonal functions for the range —1 
<K1, for 



fi 



^ _p:+,, {^)p:+,, (?)rf$= if si 5^ s, (96) 

and 



!• 



[PUsm'di- (2s+2l-\-l)T(s+iy ^^^) 

For the range 0<J<1 the set of even functions of ^, P;;+2sS, are a 
closed set of orthogonal functions, as also are the odd functions. 

If a function is developable in terms of these normal functions, its 
development is 

M)=-t:,{s+''+w J^'t^l,. PUAO (mPhs ^M^i for-l<Kl.(98) 
A special case is 



I^ 



Mi)P:+2si0d^i for 0<K1. (98^) 

For the range here needed, :<??:< 00, it is convenient to take the 
general solution of eq 87 in the form 
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where 

p..^_a + '7^) ^'^ r(M-v+l )^/ M + .'+l ^-"+1 ,3. 1 \ 

_ 2rv-H-i)i''+' 



s:W= 



Vt cos j/irrC/u+c+l) 
(l+„^)^2'',r 



Q'(i„) (99a) 



r(M-.+i)r(^ 



/I \H 2'^2' 2 2'2 '''l+W 
V2-V 



-^N^^i^^^ViW- ' (Q'-^b) 



The relation 



(i+„^)[i?;s;-i?>;i=2 cos M^rf^qr^j (loo) 



shows that these two functions are hnearly independent when 
v=^n and n=m^n. The two recurrence relations 

{2^^+i)vR;=i,^-v)R;-i-(„+,+i)R;+i (loia) 

(2,.+ 1) {l+r,')B';= - (m-^) (m+ 1)K-i- (m + "+ l)-"/'?;+i (101b) 

are also satisfied by the functions S'^. 
Since Q"{z) = iz'^—l)"D"Q^iz), the analogs of eq 92 are 

Kin) - ^~rr+;:^y ^^ (^ +^')"^^"^^(^) (102a) 

'S:(,)=^|^=^j(i+„^)"z?:s,(„). (102b) 

By use of the recurrence relations, eq 90, 91, 101a, and 101b, we 
obtain equations of the form 

(2/x+i)p;(02?;«=-AZ),[p;-i(?)i?;-iW+p;+i©5;+i(7;)] (los) 

and three others obtained by replacing RKri) by aS^(^) and then Pl(0 
by QUO ^ each of the two relations. It is evident from eq 99 



269047—41- 
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that jB^-^0 when r]-^co if ju> — 1, while SJ,~>co if /i>0. We also 
place here for reference 






and 



B'm- 



<V^+0 



-'r(^±i±-0 



(cos/xTT + COSyTr) 



(COSPCTT — COS VTt) 






andS';(0) 



(104) 



The reciprocal distance from Pii^uViA) to P{^,-n,0) has the known 
expansion 

^=xS*»cosn(e-^,) S(-l)'"(m+|) P^(0i?»('?)ntti)S::.(r;i),(105) 

provided that rj^rji, that is, the point P lies outside the ellipsoid, 
771, passing tlirough Pi, When it is inside, ??<??i and the above ex- 
pansion is merely modified by interchanging the arguments 7; and 77. 
This expansion is very similar in form to that of spherical coordinates 
given in eq 46. A comparison of eq 105 with eq 47 gives the expansion 
in spheroidal coordinates when rj^rji 



<t>n{x—xur,ri)_ 1 



V 



rri 



-^1 m=n 



S i-ir{ni+i)Pl{^)Rl{rj)Pl{^,)Sl{rji) (106) 



with interchange of 7/ and 771 in case 7;<7]^ This is similar to the 
expansion of the same function in spherical coordinates given in eq 48. 
Passing now to the case 77 = 0, which corresponds to Xi = and 
0:^ri=Ai{l — ^^)^^<iAi so that the ellipsoid, rji, shrinks into the 
circular disc, we find from eq 104 that 



SU2.+i{0)=0 and SU,M-- 



. (-1)" T{s+i) 
2« r(s+n+l)* 



Hence if 



then 



4rri 



x^+{r+riy 
<i>n{x,r,ri)_(t>n{k) 



(107) 



Vrri V 



rri 



==Z^S (^^•+^+i) r(s+t+i) ^"+"^^^ ^"-^^'^"^ Pu^{ix)> (108) 
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which is vahd for all values of 77, as is also the integral representation 

^=- ^^-l^l^J„(r0e7,(n0rf^ (109) 

-y r/'i J 

If C7is a solution of eq 64, which vanishes when x^+r'^-^co and such 
that 

is a given function of ri(=^i(l— ©|) 

on the circular annulus x={)A2<ri<^Ai, and (j==0 when o:<ri<^2, 
and when r>^i, then ^0 being defined as in eq 74 

u=2. r'?i^i%(^=2x^f r ^-%ml\/^. (110) 

JA% -yiTTi Jo ^irr\ 

Using eq 108 in this gives 

2x^, - r(.+i) 

p::+2.a). ii;:+2.(^). J%i^(?i) p'+2.($i)rf?i. dn) 

To confirm this result we note that the function TJ given by eq 111 
is a solution of eq 64, which vanishes when x^+r^ > co . 

Also Z>zC^x=±o==0 when r>^i, since {DJJ^^^={-j-B^Tj\^^ and 
P;"+2,(0)-^0. 

Whenr<^li, 27rc7--~(J9^t/)^=+o=--;^(£^,f/),=o and eq 111 gives, 
by use of eq 104, 

«^(?)--2g (2S+r.+i)j,^I^|±^ (112) 

forO<J<l. 

Keference to eq 98' for the case v^n shows that the second member 
of eq 112 is the development in normal functions of the function 
^a'(^), which is zero when Jo<f<l, that is, when 0:<r<^2, which 
shows that eq 111 is the required solution. 

The same confirmation results by using eq 109 in eq 110, which give 

L^=27r [^ e-^'''Jn{rt)di V\xcj{r,)J,{nt)dn, (113) 

Jo jAi 

which gives 

(7^-^(D,f7)x=+o=- VtJn{rt)dt V\^a{T,^Jn{ni)dT,, (114) 

^^ Jo jAi 

which is HankeFs integral representation of a function of r, which is 
equal to g when ^2<C^<C^i and is zero when 0:<r<^2 and when 
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The function 



V^l«2 



is equal to U if the density (eq 80) is used. It is then necessary to 
evaluate the integral 



r^o 1 /l-l-i^nV 2 



r(Ti+i) 

So'' 



/I — f'V 

■^r(^+l)r(l+n-A;)Jo (l-?o'cos2\«')5.sin2\J' ^^^'^■' 



Now by eq 75 and eq 81, when 17 = 



^=(l-Sf)*=(l-?o' cosV)^ so that ^,=^0 cos 1^ 



and 



5id?i=-f sin2^(^^=-^', 



SO this becomes 

Jo (l-?o' COS^ l/'jj ^ ^ ^ 

Referring to the eq 94a, we find that 

P:...(gocos ^) (-i)-2''r(s+7^+§) i-Mw^) 

(l-l^cosV)r V^r(s+l) ^ s,s+nf„,,^ocos^) 

-i ^) - ^^ ^-'*'' r(i+s-or(i+i(ra|) '^^^''' 
Now 

I cos2(i^+?i2)»/'Cos^' lA , , -±1 nis^ 

I r(^+i)rG+i) '''"2r(<+A:+7i2+-i)r(<-i-7i2+i)' ^ ^ 

Hence 



/. 



W2 

COS 2{k-^n^)^Pl+i,{^o cos i/') ^ 
(l-^^cos2^)2' 

( lV2»-> /;;y^(- IVf^' T^±s±nrH) 

(119) 
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Using this in eq 116 gives 



^num_ {-iy{\-H'.Y ^. .y,, J^{t + s+n+^) 



where 



'^'^^^^''"^"^^Sr(*+i)r(^+^+7i2+i)r(i+r~-n2-/:)r(i+7^-Z:)-^^^^^ 

In case t^ — n2, that is, i?+^2^0, 

, <-"-^>^^-M+.3+i;-(}f|)) 
'^^~ r(^-7i2+i)r(^+n2+i) 

^(~2~7 r(#-n2+i)r(^+r.2+i) ^^^^""^ 

the last form being obtained by Euler's transformation of the first 
hypergeometric function. Since t is never negative in eq 120, it is 
evident that eq 122a is apphcable whenever 722^0, and also for 
negative values of 712 when the unequality C>—n2 holds. In the 
remaining case, where 0^t<i—n2^ we may write ?i2— — m, where m 
is a positive integer, so that, 0:^^:^m. 
In this case 

=r(«,+ i)(-i) \p^j r(2<+i)r(«+7^-TO+i)f(^-i!+i) 

(122b) 

This expression however is that which eql22a becomes when i+W2<C0, 
so that either eq 122a or 122b may be used, as they are both valid 
for positive and negative values of n^. Using the latter, where m is 
replaced by — «2 gives 



(-1)V5',== 

( 



i)"{^?) r(.+ i) ,^,^.;-^.^^^^,)r(2^+i)r(i-..-o - ^^''') 
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Using this, eq 120 becomes 

s 2^r{t+s+n+^)F(^~2t,n+l,l-n2-t;'^'^ 

^r(i+s-or(«+w+7i,2+i)r(2i+i)r(i— Mo-O 

^(-l)'+ "Hl + ro)"+"'(l-g ;)-"'^^ 2''r{t+s+n+i) 

2t'^^AI HT{l+s-t)r{t+n+n2+l) 

^T[k+l}T{k-t-n2+l)T{l+2t-k) ' 
Interchanging the order of summation gives 

2^5^^ ^^ ^\ 2/ r(A: + i) ^*' 

(124) 
where 



<=o 



r(2H-i-^)r(^+^+r^2+i)r(i+s~0r(i+A:-7^2-0 



In order to evaluate Gjc (eq 125) as a single term, we define the 
function G (z) as follows: 

^ A A {-2y2'*T(t+s+n+^) 

kt:^o,^Ti2t+l-k)T(t+n+n2+l)Til+s-t)T{l + k-n,-t) 

(126) 

Since n and 712 are integers, this function is a polynomial in positive 
and negative powers of z. It is therefore analytic at every finite 
point in the s-plane except at s=0. After inverting the order of 
summation, we may substitute k=fi+n2+t and sum first over ^ 



7 W ^ (~z)''^+'2"rit+s+n+i} 

^^' ^T{t+n+n,+ l)ril+s-t)rit-n2+iy 

- {-zyr{t~n,+i) 



^f='„r(f-n2+i-M)r(M+i) 

or 



. ' 2^'r(f+s+Ti-+-|)(-g)"'+'(i-g)'-"' 
^^^""^or(«+n+7i2+i)r(n-s-<)r(<-n2+i) 



(127) 



(128; 



= 

This series starts not below t=n2. Let X=<— m2, then 



^f.-\ ( o-^2.^v^ [-42(l-g)?r(X+s+Ti+n2+i) ,,2q^ 

^(.) = (-2.) 'gr(x+i)r(x+7t+2n2+i)r(i+s-n2-x)- (12^) 
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The sum is of the form 

By a special transformation due to Gauss 
and by use of the equation 

_ 2ra-^) ^ (-g)^r(X+2a) 

We get 

*= — CD 

_ (-2)2"2r(s-7i2+i) ^ (-2)^r(X+2s+2n+2w2+l) 



22»+-''2r(s+rj,+ri.2+l) iSo r(X+l)r(X+n+27i2+l)r(l+2s-27!,2-X) 

(131) 

from whicli, by equating equal powers of (z), we get for G/c with k = 
X+2712 the expression (eq 132) below. This brief method of deriving 
eq 132 from eq 125 was suggested to me by C. L. Critchfield, whose 
helpful criticism in other respects may be acknowledged here. 

P.7. , „ , .^^ T(s^n2+i)r(k+2s+2n+l) 

Using this in eq 124 gives 

rn,.n. (-i)^+^^r(^-7i2+|) / 1 + ^;\ Vi " g; V ^» 

^?7r'/»2^+2'^2+ir(5 + 7^+n2+l)\ 2 /•yi + U 



2s /l—t'Xfc 



r(^+2s+27i.+ l) 



T{k+l)r{k + l-2v^)T(l+2s+k) 
Now 

^r ,^/l-?iY T(k+2s+2n+l) 



§<-"*(-/) 



Tik+l)r{k+l-2'n2)T{l+2s-k) 
T(2s+2n+l) 



"r(2s + l)r(l-2n.)K-2«'2« + 2.+ l,l-2^;y») 

Kj;) ^('/yimill) ^^ ^^33) 

by Euler's transformation, remembering that 1— ?o==2a2Mi and 
l+$o=2a,Mi. 
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We thus obtain the following 

where the ^-function is defined by 

_ r(s+rii4-i)r(g+7i2+i)r(s+n.+TO2+i)r(s+i) 

~ r(2s+l)r(2s+2ni + l)r(2«,2+l) 

On expressing this as a polynomial and summing it in the reverse 
order, we find, since s is an integer, 

AlT''''Kai,a2)==AlT^'''^(a2,a,). (136) 

Also since 7^2 is an integer, it is found that 

A.^"'-'"'(an,a,):=A^i-lZ(aua,), ' (137) 

which vanishes if s—n2<^0. 

If n2=ni in this, it is found that 

\2ni~2ni(^ ^\_ K-2nu2ni(^ \_. ^F (^ + |) T (s — 7^1 + j) _ip2ni(^'\ (^oq\ 

^28 l^b"'2j— A2s ^ ^^ '^^~2^'^^V{s+l)V(s+n +1) ^ ^^ \^^^) 

which vanishes when 5<?ii. 

Recalling that 7i=ni— 7i2^0 and ni>:0, it is evident that eq 134 may 
be written 

In this equation 0^'''^^ is the value of hPn-\-2s(^i)<^i^i)d^ii 

where the expression (eq 80) has been used for the density o", so that 
with this value of the integral the formula (111) is the expansion of 

the function i!7=<^ni,n2fe^jCti,a2)/V^i^2. 

We have, therefore, obtained the following development when 
n=ni—n2^0j which is valid for all positions of the two parallel circles, 
which converges to a finite value except in the case where the circles 
are equal in radii, and approach coincidence (in which case it becomes 
logarithmically infinite). 

0npn2(^/>^l;<^2) 

==^^^i2(-^yi^s+n+i)Al^^'-'-\a,^^^^^ (140) 

7«'-^l 8 = 

When n9.=%^0, this becomes by use of eq 138 
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<t>ni,ni{Xyr,ai,a2) _ 

From eq 137 it is seen that the expansion (eq:<140) begins with the 
term s=n2 when 0:^n2^nij but with s=0 when rZi^O and ?Z2<0. 

It may be noted that if n^\ni—n2\ and if we write A^^'^''^^^''{ai,a2) 
in eq 140, the expansion is vahd for all integral values of % and r^2, 
with the understanding that the upper of the double signs is to be 
taken when ni—n2^0, the lower when %— '//.2<0. The ambiguity 
disappears when 7ii — n2=0. 

From the integral representation (eq 57a) of ~j=LJ. it is evident tliat 

-y/aiao 

this function is unaltered if we interchange the arguments ai and a2 
and corresponding parameters ni and n2 and multiply by (~l)"i~'^2^ 
This fact is derivable from eq 140 by use of the relation eq 136, 
noting that Ai, A2, f, and rj are all unalfected by interchange of ai 
and a2. 

The integral rej^resentation shows also that 

4>n^',n.i{x,r,ai,a2) __(—!) ''^4>n^yn2{x,ai,r,a,^) 

V<^1^2 V^'^2 

Where 



^1 = '%— ^2 



(142) 



This interchange of arguments r and aj in the light of eq 65 amounts 
to viewing x and ai as the cylindrical coordinates of a point, while r 
and a2 are constants playing the part which ^i and 02 held in the 
preceding formulas. 

Hence, as in eq 173, we let 

A=r+a2, x=A,1^, ai-liVa-?)(l+?), (143) 

so that f and rj are computed from the second members of eq 82, after 
replacing Ai by r+^2 and r by ai. 

The function on the right of eq 142 is then given by a development 
obtainable from eq 140 by replacing Ai, ai, ^, 77, ni, 7^2, and n by Ai 
^i 4*j 'Hi ni—n2j —n2, and \ni\, respectively. This gives 

0;,j,„/x,r,ai,a2)_ l 



Vai^s 7r''"(/'+a2) 

±{~m'28 + \n,\+^^)AI^^^-.-n,^^^^^^^ (144) 

the upper or lower sign being taken according as ni is positive or 
negative. Similarly, from interchange of r and a2, 

<i>nim2^Xyr,ai,a2) ^ <j)nji^(x,a2,ai,r) 

'y/aia2 -y/aiT 
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where 

(145) 



n^^nx—ni 
n^\ni~n'^=\n'^ 

This gives a third expansion 



<l>ni,n2(^,r,ai,a2) _ (—1)711-712 

V^ ~7r^''2(r+ai) 

g- v(^2s+\n,\+iy2s(r,a,)'^''^^^^^^^ -RIuUsQ^ (146) 

where the upper or lower sign is taken according as ris is positive or 
negative 



x={r+a{)^rj anda2-(r+ai).V(l-«')(l + ^'), (147) 

so that I and r) are computed from the second members of eq 82 after 
replacing Ai by r+«i and r by a2. 

The following section is concerned with the partial integrals with 
respect to x derivable from the expansion (eq 140). 

VII. EXPANSION OF c^n.m IN SPHEROIDAL COORDINATES 

The function coni,n2(^/,«i,«2), or more briefly, oini,n2(x'^) is represented 
in eq 59 by the integral 

co«,.w2(^^0=4^V^i^2 dx'' (t>n,,n2(x',r)dx\ (148) 

Jo Jo 

We exclude for the present the exceptional case where ni and 712 are 
both zero. In all other cases we may write 

oini.nii^yr) = Vn,,n2{x,r) — Ki.n2(0/) +27r|x|S';,j.n2(/') (149) 

^'n,,n,iXy'd==ynun,{x,'f')±2TrSnun2{r), (150) 

the upper or lower sign being taken according as x is positive or 
negative. The function F'«i.«2(2^/) is an even function of x defined by 

T^«^n2(^/)=47^2ala2J^6-I^^'J;,l_^2(r0«/;^l(al0e7;^2(^20^• (151) 

The function V'ni,ni{xir) is an odd function of x defined by 
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the upper or lowei 
negative. Also 


sign being taken according as a 


; is 


positive or 


STrSni.Tpg W = "~ Vni.n2 


(+0,r) = + F;,„,(-0,r) 










=4ir^aia2 1 Jni~n2{rt)J„^{ait)J„^ 


1 *\^* 


(153) 


These relations are all obtained by introducing into 
tegral representation (eq 57a) for (I>nin2j which gives 

r 


eq 


148 


the in- 



If the exceptional case 7ii=7i2=0 be excluded, this may be separated 
into three integrals, each of which converges, thus giving the eq 149 
to 153. 

From 152 arid 57a it is evident that 

1^711,712(^,0 = —4:Tr^-y/aia2 I (t>ni,n2(:^',r)dx' (155) 

if x>0, 
and 

if x>0. In both of these integrals, the integration must be performed, 
while r is held constant. 

The functions Vni,n2 ^^^ ^7ii,7n ^^^^ satisfy the partial differential 
eq 64. The function y„',,„j vanishes when the point P(x,r) goes to 
infinity in any direction in the (x,r) plane, while Vn^,n, vanishes in 
general if x-^dsz «» , while r has any fixed finite value. 

To perform the integrations in eq 155 we write eq 103 in the two 
forms : 

-~^PJP:+2.-iI?) •S:+2.-i(^) + H+2.+i(?) •i?n+2.+i(^ (157) 
(2s+n+f)P^+2.+i(^)^"+2.+i(r;) = 

-^D,K+2.(s^) • K+2sM+P:^2s+2(^) ' R:+2s-^2(v)1 (158) 

The expansion for Vi^^m i^ obtained by using the expansion (eq 140) 
for (j)n,n2 ill <^q 155 and then performing the integration with respect 
to X by use of eq 157. Since r is held at a finite constant value in 
this integration, the value of the integrated terms at 3:== 00 (r finite) 
are found by letting ry-^oo, while at the same time |^->1, the two 
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variables being connected, as shown in eq 81, in such a manner that 
Umit 

It is thus found that for the case n=ni— n2>0 



+s(-i)'WL"""'"'(«i.<^2)-^?r^^'"'(a.,a2)]p:;+2,+i(i)-B:+2,+i('?) (159) 

5=0 J 

and for the case n2=n{^0{n=0) 

V' Ax, r, ai, a^) = (- 1) ^^^'2a,a,^T ± (- 1)^^^/"''-^"^^!, ^2) 

" 5 = Til — 1 

-^2^+2""'"'"' (*" (h)]P2s+liQR2.+i{v). (160) 

The first term inside the brace of eq 159 is 

=^^i^^¥i^.f«i^(l,.+i,H-l;l-a, (161) 
for from eq 95b we find 

r P::+2,-i(Q 1 __21p(^_j_i)(3_j2)-2.^2^(l^ ,,+^^ ^„f_l. i__^2)^ (i(32a,) 
L ^ Js=o n-^-K 

and from the eq 99a we find 

n 

W.^,.-M)\.-.-^f^^- (162b) 

Also from eq 95a we find 

Gp;.»-,tt)}..=- ^-'""':|"-"' ^(i.»-i.'.+i;'-«. 

and from eq 99a 
Hence 

[p"+2.-2(i)..B:;+2,-2('?)L=o=^^i^2y^(i,"-i.^+ 1 ;i -?') • 
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Inserting the expansion (eq 159) in (eq 156) and integrating by use of 
(eq 158) gives for the case ti^O. 

I nin—^j-ylTT \l + 77 / 

^--^^^^^^^fq^^^f^^^ (163) 

and for the case 7i=0, ?i2=ri-i>0 



9q ] 



The value of Vnx/a^ {0,r^ai,a^ for the range 0^/':<^i is found by 
placing ?7 = and J=^/i_Z-in eq 163. For the range 7*>^i we 

place ^=0 and 77=^ /^2~~1- Similarly, from eq 159 the value of 

2TSnirii(j')~ — Vn,n2(-\-0,rjaia2) could be found, but such expansions 
are unnecessary, for we now find that Smn^i'^") may be obtained in 
finite form. 

It should be noticed that the term ^t\l,n+l/2,n+l] 1 — ^^) in 
eq 159 is an odd function of ^ (and therefore of x), as are all the terms 
in Pn+2s-\-i{0 in the series. The latter all vanish with ^, but the former 
does not, as may be seen by the transformation 

iF{l,n-\-hn+l;l~-e)^±F{hn,n+l; 1-^- ^^^J_i'\i'^^- (165) 
when ^->it:0. 

Hence when x-> + and r>u4i, i. e., J-->-l-0 and ^ — ■ — r,^-S eq 159 
gives for r>^i 

if n^^O. Referring to eq 135, we find 

V-r(n+i) ^Ar(i+^i)Ar(i~7^2)/ ^'""^^ 

when ni>0 and 7i2<rii. 
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Hence eq 166 becomes 

when r^Ai=ai-\-a2 and ?ii>0, n2<^ni. This vanishes if n2 is a positive 
integer less than n^. 

This relation may be used to obtain the value of Sn^rn^f) for the range 
Q-^T<^A2 by considering the effect of interchanging r with ai or ^2, 
which is determined by the integral definition (eq 153 of Sn.mi^)^ 
which is an abbreviation for S'^in3(^,<^i,^2). It is thus found for the 
case 7ii>0 and n2<% that when 0:<r:<A=|ai— ^2! 



S'^i'«^(r)=:.7raia2: 
= 



r(ri+l)aTr(n2+l) 



, when a2<Cai 



when a2>(Xi 



(169) 



Similarly, for the case n2=%>0 it is found that 



Sni,ni{r)=^0 when ^i<r<oo 
ni \aij '^ 



when 0:<r<A. 



(170) 



In all these statements the case ni=n2=0 is excluded. It is evident 
therefore that Sn^^mi'f') vanishes when /•— 00 and also iini9^n2 when 

ni,«a v) IS a 

To find 



r = 0. It may be shown from its definition (eq 153) that S, 

continuous function of r, in particul^ at r=Ai and r=A2. 

its value for the range A2<Cr<^Ai, we find from eq 152 that when 

x>0 

t—D (r'T)\ '^''^ Y^nx,n2\Xir) YX2 ^nun2\Xjr) 
^Ur\riJr) ^2j 27r ^^ 27r 

= 27raia2 I te-''^Jn^_n2iri)Jn^i(^\t)Jn2(S^2t)dt. 

Letting .t approach zero in this, and noting that by eq 153 
y'nu'nA^^^'f) ^ — 27rS'ni,n3(/'), and that the integral on the right approaches 
that which defines 27ro-(r) in eq 66, we find that Snumir) satisfies the 
ordinary nonhomogeneous differential equation 



1 



Dr(rDrSm,m) ' 



n'S., 



^l'«2 _ 



-47r^aia2(7(r) 






cos (tii^?— 712^2) 
sin (6>^-6>;) 



(171) 



for the range ^2<C^<C^i> this value of a being taken from eq 72. 
If r lies outside this range o-^j^^^O, the homogeneous equation having 
solution Cir'^+Cor''^ ii n±0 and C1 + C2 log r if n=0. The values just 
found are of this form. Hence if any particular integral of eq 171 
is found for the range A2<ir<^Ai, this increased by Cir"" -\- C2r~ "^ is the 
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general solution for that range, and the values of the two constants of 
integration, C^ and 62, are then determined by the fact that this 
solution must be continuous at r=A2 and at r=Ai with the solutions 
previously found and given by eq 168 and 169. Similarly, in the 
case 71=0. 

To derive a particular integral of eq 171 in finite form for the case 
7ii>0 and n2<Cni, so that 7^=ni— 7i2>0, we begin with the elementary 
integral form of the solution 



S{r)=^[r-''I,-r-I,], (172) 



where 



/i=- (r'^cos {nie\—n2el)d2\l/ 

I,= j r-'^cos (nidi — n2dl)d2\P 

Now by eq 77 

r^ cos (r^i^;~7?.2^^)=^<6-2^^^2/i_?2^\ 

^(-l)^(lJ cos 2(^+7.2)^ 



2<^\ n 



=arr(?i+i)S 



^0 Tik+l)T(l+n-k) 



= (-i)"24r(n.+ i)Z; 



{-iy(^y cos2krl^ 



so that 



a^' ' 'i^,T(,k-n2+l)T{l+n,-k) 



r A-( 1) ^„ a, a, '^(r(n,+i)r(i-7i2) + 



<-»'(S)' 



sin 2kyl\, (173) 



^n,V{k-n.+ l)V{l+n,-k) k J 

where S' indicates that the term Ar=0 is omitted. To evaluate I 2 
we derive from eq 77 



and 






so that 



ai 
where the real part is to be taken 



■Sw'/(?)[' +?]'<?> 



r re^Tt'- 
Since ni>0, the expansion of the binomial \-\ is a finite 
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polynomial in re»°a2. This gives 

_r'T(ni)a2"2r go 



r"L 



ai" 



r(n)r(ft2+l) 



m 



n-l 



sill k6l\, 



(174) 



Using eq 173 and 174 in eq 172 gives a particular integral S(r) of 
eq 171 for the range ^42<r<4i, when ?Zi>0 and ?i2<ni. Hence if 
A2^r<CAu %>0 and n2<7ii, we find 



= ai% 



(-i)^r 



r- [T{n, + 1)T{ 



)r (1-712) 






(-1 



'■(ST 



<■ Lr(n+] 



/t=n,r(/t-n2+i)r(i+fti-i-) 
el 



sin 2^1/- 



l)r(n2+l) 



1 ^""^ 



sin A:i9S 



(175) 



where the constants Ci and C2 have been so chosen as to make Snimir) 
continuous at r=Ai and at r=A2, its values at these points being 
given by eq 168 and 169 when r=Ai or A2, respectively. 

The terms in the summation vanish at each of these values of r, 
for, when r=Ai, 61=0, ^2^7r, and 2ip=ei—e[==7r, and when r=-^2, 



2iA^ 



^0, andf;hJh^if-^^<^^ 



This shows that the second member of eq 175 reduces to eq 168 
when r=Ai and to eq 169 when r=A2. Therefore, eq 168, 169, and 
175 give Srnin2 in finite terms for all values of r when %>0 and 7i2<C'^'i- 

In eq 175 Sn^n2(T) is Sn^n2('^,(^i(^2)' In order to find Snin2('^,a2,ai), 
we must not only interchange ai and ^2 where they occur explicitly 
in the second member of eq 172 but must also replace 62 by tt—Oi, as 
shown by eq 78b. 
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The remaining irreducible cases are for n2=ni^0 and n2=ni=^0. 
In the case ng— ni>0(7i^0), we start with the particular integral 
of eq 171 

S(r) = 4:n%a2 j no-(n) iog^dri=2aia2 I cos 2ni\p' log -6^2^' 

=aia2^ ^ sm 27^1^ log ^+2 fcos 2n^r log nd2r\ 

= a,a2[^^f^log^^-2±j^^^^ 2n.^^ cos 2krd2r] 

^a4?-?^^^^log^^ 
I Til ^ r 

This may be reduced to finite terms by use of the expansion 
2 sin 2n^\p , ?1^^}Y^2Y2 sin 2ni\p cos 2^1/^ 



ni ^ r fc^A:\ai / Ui 



so that 

S(r) = ^'^^^y^f^AT ^^^ 2 (^ + Ui) \l/ __ sin 2(lc—ni)\p 'l 
nii^AaiJ[_ k+ni k—n^ J 

-"-fl-0"'-[(:i)"-(s)"'H:-;)'^ 

yKHMV"''^ /aiV'^'n sin 2k^p sin 2uif 
Now 

^ie^^l ^ 

\ Gil 
so that taking the logarithm of both sides we obtain 

e?=g(^'J'-^^- (177) 

The particular integral then becomes for the case ?i2=r?i>0. 

^ ^ ni 

269047—41 4 
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which contains only a finite number of terms and is vahd whichever 
of the two circles is the larger. 

With this particular integral we find, after determining the two 
constants of integration by making Smm ir) continuous at r=A2 
and at r=Ai with the values given in eq 170, 



&.«. W ^^ y^^ 



1 — ) II a2<Cai 

ni \aj '^ 

= (— ) ifai<a2 

-°fi(i)> +©"■('-«■ 



when 0:^r:<^2? 



sin' 27^1 1/^ 

Til 



sin 2kyp ] 

when ^2- 
when Ai- 



lA, 



= 

provided that 7?.i>0. 

In the particular case of greatest importance, ni=l, 

\ Tral if a2<<^i 
Wal if ai<^a2_ 

= ale i+aKir— 62)— aia2 sin 2\l/ when A2^rz 

=0 when Ai^^rz 



(178) 



Sn{r)- 



:r<^2 
:^1 



(179) 



This shows that Sn(r) is the area common to the two coplanar 
circles of radii ai and a2, when r is the distance between their axes. 

By the preceding formulas every o)n^n2 is evaluated with the excep- 
tion of the case ni=n2 = 0. Before considering it, attention may be 
called to the interesting fact that although every T^wpna vanishes when 
x->± 00 holding r finite, there are two cases of eq 163 namely, 
Vi,o(Xjr,ai,a2) and Fo^-i (x,r,ai,a2) which do not vanish when the 
point P(x,r) goes to infinity in directions other than the a:-direction. 

By use of eq 163 and 167 it is found that when x^+r^-^oo in the 

X 

direction (j) defined by cos <^= , ^ . == then 



Vuo(x,r,ai,a2) -^ 27r^ala2 
Vo,-i(x,r,ai,a2)-^—2Tr%al 



tan 



<!> 



tan I 



(180) 



In the exceptional case where 712= ^i = 0, the integrals involving 
BessePs functions, which served to define the functions Vnum and 
Vjium in eq 151 and 152, are not convergent. Neither are the defi- 
nitions of these functions by the integrals, eq 155 and 156. Neverthe- 
less, functions 

Voo{Xjr)j V!io{Xjr) ^Dj^Vooix^r), and Sooir) 
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may be defined in this case, so that eq 149 and 150 are apphcable. 
To find these functions we begin with 

cooo(^/)=- {\oo{x',r)dx'=4.Tr'a,a2 i^^dx' P h^^;^ dx' (181) 

Jo Jo Jo V<^1^2 

coiofer) =D,c^oo{x,r) =47r'a^a2 ('^^d^ dx\ (182) 

Jo V^1^2 

Applying the operator r~Wr{rDr) to the latter gives 
r Jo r \ yaiaa / 

Jo \ V^1^2 / 

\ V^1^2 / L \ V«1^2 / J^->0 

Also applying the operator Dl to eq 182 gives 

V^1^2 

so that 

vWo(a:,r) =fe+iz>.(ri5.)1a;o'oC^/) =:4^^aiarO,%M1 

That is, by eq 66 

V^Woo(:r,r) = =F 4x2^1^2. 27ra-oo(r) (183) 

the upper sign if a:>0 the lower if a:<0. This shows that a>^(a:,r) is 
of the form as in eq 150 

coo'o(:r,r) = Voo{x,r)±27rSoo{r), (184) 

where VQo(x,r) is an odd function of x satisfying 

V-Fo'o(a:,r)=0, (185) 

and Soo(r) is a solution of 

V'^oo(^) =lDr{rDrSoo{r)) = -4.T\,a,aoo(r) =^^> (186) 

as ineq 171. 

Since cojo and coqo both vanish with x^ 

S^ir) = -Y^^^^+Y^^^, (187) 

as in eq 153. 
Also by eq 184 

o^oo{x,r) = j^ccoo{x',r)dx'==Voo{x,r)-Voo(0,r)+2w\x\^^^^ (188) 

as in eq 149. 
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If in the integral (eq 182) we introduce for 0oo(^^^)/V<^l^2 its ex- 
pansion (eq 140), where -^2^=^27^^ P2s{Q by eq 138, this g: 



ives 



w'oofeO = 



^^ala2^h 



00 K-^)' J A 



where 






(190) 



The terms in the series of eq 189, for which s>0, may be integrated 
by the use of formula 157, but the present case is again exceptional 
for eq 157 now fails entirely, for the term s=0 becomes infinite instead 
of merely becoming indeterminate of the form 0/0, as in the case of P;', 
where n=0. 

This breakdown of the formula may be repaired by writing 



dfx 



Pott)i?o(.)-27j.=o^^^^)^^^^)7 



(191) 



where the integral is a contour integral taken around a small closed 
curve in the complex ^i -plane, encircling the point ^==0. Since ju does 
not vanish at any point on this curve, the integrand may be written 
by eq 103 

which gives by Cauchy's theorem 

Po{mM=j^ cot-' r,= -AJ)\D ^^^-^f^^'' ^''^ 1-0+ ^' (^)^i(^)}- 

(192) 

(193) 



It is found that 



V 



;rr, V^P,-x{m.-i{-n)l i^„. (l±|)VI±r. 
r^l 2m+1 J ^ 2~ 



If we define the function v{x,r) by 



«(»;,/•) = ^log 

= + 7= log 

sjTt 



(! + ?)• Vl + ^^ 



{l-&^l + ri' 



when a;>0 
when a:<0 



(194) 



then v{x^r) is an odd function of x, which is an axially symmetrical 
solution of Laplace's equation, and eq 192 becomes 

P,{mM-=^<^oi\^-A,BMx.r) + P^{mi{'n)]- (195) 
A further integration may be made by the same method. This gives 



2A 



(196) 
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The function within the parentheses of eq 196 is an even function of a;, 
satisfying Laplace's equation. Integrating cq 189 by use of eq 157 
and 195 leads to the following definition of ^00 

V^{x,r)=^—2aia2-y/Tr 

[7iMx/)+^^i~mA.Ts{aia2)-A^Uai^^^^^ (197) 

Integrating this from x=0 to x=x by eq 158 and 196 leads to 

Voo(x,r) ==-~a,a2^\2Axv(x,r) +^l v]+A,J2(- ly 

L 2^"i 2H=1 J^2,(^)i^2.Wj- (198) 

From eq 197, 194, and 187 we find 

So(i{r) = — 27raia2 log^-^ When r>^i. (199) 

When 0<7'<A, Sooi^) is found from eq 187 and 197 by placing 

/ P" 
77==0 and ^=^^-=^1—21 

Soo(r)-^-^[47rlog(l+?.) 

+;^S (- 1)^ i^^l [^2^-22 (a,,a2)-'AT.(an,a2)]P2s,.i (^)) (200) 

for the range 0<r<^i. In the part of this range 0<r<yl2 where 
o-QQr=o, this development must represent 

Soo(r)=Co+Ci\ogr, (201) 

as shown by eq 186. 

If the development eq 200 converges when r— 0(^r=l), the constant 
Ci must be zero. 

Now by eq 200, when r=0, Jr-=1 and P2s+i(l)==l. Also reference 
to eq 138 shows that 

Aoo ttF (gH--^-) p /«./\ 

-^2s r2/'^ I 2^ ^2s\Koj» 

Hence from eq 200 we find 

Soo(0)--27raia2{log 2-F(Q}, 
where 



nQ-i+j-^± (-i)»^[i+2^] A.(io: 



;). 



Since this series converges, the constant Ci in eq 201 is zero. It is 
readily found that F(^o) is the development 



««) = -%[i^]=-log^ 



-} 
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so that 

Sqo(0) = — 27raia2 log ^^ ^ where —^ — ^=ai if (i2<^i 

= a2 if a2^ai. 
Hence eq 201 becomes 



when O^r^A. (202a) 



Smif) = —2iraia2 log ( ^~ ) if «2<c^i 
^-27ra,a2log(^^Jifa2>a, 
Also for the range A2<Cr<CAi, we have 

Soo(r) = 2a.a.jW log ^+. log (^)-g(|y ^{ if a.<a. 

^2a,a2\2^ log -+. log (^_--- j-g(^- j -^^j if a2>a. 

(202b) 
Finally, for the remaining range of r, 

Soo{r) = —27raia2 log ( , j when ^i<r<oo . (202c) 

This completes the expansions required. Before summarizing the 
results, we should consider a small correction which takes account of 
the diameter of the wires. 

VIII. CORRECTIONS FOR THE DIAMETER OF THE WIRES 

If a helix is made with ordinary wire, a plane section of the wire 
by a plane thi-ough the a^-axis may be considered circular. Let pi 
and p2 be the wire radii, and let ai and a2 represent the cyHndrical 
radii of the central helical filament of the actual helices. Since pi 
and p2 cannot be larger than irpi or Tp2 the correction terms for finite 
thickness of the wire will be second-order infinitesimals. 

Assume that the current density at any point in the wire is in the 
direction of the generating helical filament passing through that 
point, its magnitude being a function of the distance r, of the point 
from the axis of the cylinder on which the wire is wound. Let ji 
(or 72) be the value, at the central filament, of the r-derivative of 
the current density. 

It may then be shown that the correction for the thickness of the 
wires may be made by using for the cylindrical radii ai, a2 in the 
principal term not (Zi, 02 of the central filaments of the wires, but 
ai, a2 of their equivalent helices, where 

ai = ai+|(^ + 27i) and a2==a2+|(^^+272)- (203) 

If the current is uniform, 71=72=0. If it is the ^'natural" distribu- 
tion in both wires (inversely proportional to r) then 71= — IM and 
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72= — 1/<X2, SO that the correction is equal in magnitude but opposite 
in sign to that of the uniform distribution. 

The derivation of eq 203 is almost verbatim, the same as that 
given in the paper on coaxial helices (footnote 2). It depends essen- 
tially on the fact that ^n satisfies the partial differential equation 65. 

IX. SUMMARY 

The mutual inductance, M, between two parallel helices is ex- 
pressed in eq 27 in terms of a function w (x, 61,62) involving four sets of 
terminal ^-differences defined in eq 22 and 22', where 61 and 60 are the 
associated terminal azimuths, as given in eq 29. In eq 33 and 58 
£0 is resolved into four others, ccii-{-pip20)oQ/aia2+o)i + oi2, the first 
two, coii+2^i^2WooMi^2, corresponding to the mutual inductance of the 
two associated current sheets, the principal term con representing the 
effect of angular components of current, the relatively small term, 
jPiP2coooMiG^2, of second order representing their axial currents. The 
terms coi and C02 depend upon x and the terminal azimuths. These 
are expanded, eq 61 and 62, in Fourier's series of the terminal azimuths, 
the coefficients involving functions o)\^,n2 and <^nin2- It is shown 
that two terms of the series for the first-order quantity coi are sufficient 
in general, the second-order term C02 being given by retaining terms 
in the single and double series corresponding to ^1 — 1, 2 and n2=ly 2. 

In eq 140, 144, and 146 three different types of expansions for any 
function ^^ng ^^ oblate spheroidal harmonics are given. In section 
VII the first and second partial integrals of these expansions are 
evaluated, giving the functions o)\i^n2 and o)ni,n2 ii^ spheroidal har- 
monics. The results are found in eq 149 to 202. Finally, in eq 203, 
is given the effect of diameter of the wires. 

Washington, August 8, 1940. 



